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1 Evidence for a hemispheric asymmetry in basin sizes 

In our analysis we considered large impact basins that possess unambiguous signatures of 
crustal thinning in the GRAIL-derived crustal thickness model (9). Basins with crustal thin- 
ning diameters, D, less than 200 km were not analyzed in order to minimize preservation biases 
that might result from burial by mare basalts or ejecta from adjacent basins, or from deficiencies 
in the crustal thickness model. Azimuthally averaged crustal thickness profiles were generated, 
and the diameters of crustal thinning determined from the radial distance from basin center at 
which the crustal thickness first reaches its pre-impact value ( 1 , 2) (Table SI). The pre-impact 
crustal thickness, H, was taken as the average value at a distance of three basin radii from the 
basin center. 


Table SI. Location, diameter of crustal thinning D, and pre-impact crustal 
thickness H, for basins with diameter greater than 200 km. Stratigraphic ages 
correspond to the Imbrian (I), Nectarian (N), and pre-Nectarian (PN). 


Basin name 

(°E, °N) 

D, km 

H, km 

Age (3) 

Relative age (55, 36) 

Nearside 

Imbrium 

(341,37) 

680 

28 

I 

3 

Serenitiatis 

(19,25) 

607 

28 

PN 

4 (uncertain) 

Crisium 

(59,17) 

476 

31 

N 

4 

Smythii 

(87,-1) 

465 

35 

PN 

11 

Nectaris 

(34,-16) 

434 

34 

N 

6 

Humorum 

(321,-24) 

380 

31 

N 

4 

Asperitatis 

(27,-8) 

370 

34 

- 

- 

Humboldtianum 

(83,58) 

325 

31 

N 

4 

Cmger-Sirsalis 

(293,-16) 

270 

35 

PN 

1-7 (uncertain) 

Humboldt 

(82,-27) 

235 

31 

- 

- 

C25 (57) 

(351,11) 

235 

34 

- 

- 

Schiller-Zucchius 

(315,-56) 

209 

35 

PN 

9 

Farside 

Orientale 

(266,-19) 

418 

43 

I 

1 

Moscoviense 

(148,26) 

319 

42 

N 

6 

Freundlich-Sharanov 

(175,19) 

308 

44 

PN/N 

8 

Coulomb-Sarton 

(237, 54) 

305 

42 

PN 

11 

Fitzgerald- J ackson 

(192,25) 

303 

49 

PN 

- 

Mendel-Rydberg 

(265,-50) 

286 

41 

N/PN 

6 

T22 (57) 

(179,49) 

281 

45 

- 

- 

Apollo 

(209,-36) 

276 

24 

PN/N 

7 

Hertzsprung 

(231,2) 

252 

54 

N/PN 

4 

Compton-Belkovitch 

(105,61) 

230 

29 

- 

- 

Korolev 

(202,-5) 

223 

57 

N/PN 

6 

Dirichlet- J ackson 

(201,13) 

214 

56 

PN 

- 
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The cumulative number of basins with diameters greater than D is plotted in Fig. 2 for both 
the nearside and farside hemispheres. We note that there are 12 basins with diameters greater 
than 200 km on each hemisphere, and nine basins with diameters greater than 250 km on each 
hemisphere, consistent with the hypothesis of a globally uniform impact flux. Nevertheless, 
Fig. 2 shows that the size distributions for the two hemispheres differ significantly. With in- 
creasing diameter, the number of basins on the farside drops considerably faster than that for 
the nearside. For D > 320 km there are eight basins on the nearside and only one on the farside. 
The probability that such a distribution of basins would occur by chance under the assumption 
of a uniform impact flux is quantified in Fig. SI. The total number of basins with diameters 
greater than a specific value was first determined, and the probability that the number of basins 
on the nearside would be greater than or equal to the observed value was calculated from the 
binomial distribution. For diameters greater than or equal to 320 km, there is less than a 2% 
probability that such a distribution could have occurred by chance. Though the probability in- 
creases at large D, this result is simply a reflection of the fact that the total number of basins 
decreases with increasing diameter. Previous studies of the Moon’s impact bombardment by 
near-Earth asteroids show that the cratering rate between the nearside and farside hemispheres 
should be similar {12, 33, 34) for a large range of impact conditions. 



Diameter (km) 

Fig. SI. Probability that the number of basins on the nearside hemisphere of the Moon with diameter 
greater than D would be as great or greater than the observed value. 

The impact basins with crustal thickness signatures in this study correspond to approxi- 
mately the youngest half of the entire inventory of lunar basins. Among the 15 relative age 
groups of Wilhelms (35) these basins are all younger than or contemporaneous with group 11. 
Globally, there is a total of about 27 basins that formed with ages less than or equal to group 1 1, 
and about 19 that are older {36). Stratigraphically, seven of the basins with crustal thickness sig- 
natures formed in the Imbrian and Nectarian periods, five during the Nectarian or pre-Nectarian, 
and seven during the pre-Nectarian (3). The youngest basin in our study, Orientale, is estimated 
to have formed 3.73 Ga (billion years ago) {12) and an upper limit on the beginning of the 
Nectarian period is about 4.2 Ga {11). The South Pole-Aitken basin was excluded from con- 
sideration in our study given that it is the oldest impact structure on the Moon (it is the sole 
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member of group 15 of Wilhelms), and because the thermal regime of the crust and underlying 
mantle at the time it formed is highly uncertain. 

In addition to the impact basins utilized in this study, we note that there are three large 
basins with ambiguous crustal thickness signatures: Australe (94° E, 34° S, 17=6003=300 km), 
Nubium (343° E, 21° S, Z3=700±200 km), and Fecunditatis (52° E, 5° S, 79=450±250 km). 
These three basins are all located on the nearside hemisphere and are older than those basins 
with known relative ages in Table SI. Inclusion of these basins in our study would have made 
the asymmetry in impact basin sizes even more pronounced. 


2 Impact modelling setup 

iSALE-2D (21-23) is a multi-material and multi-rheology finite difference shock physics code 
used for simulating impact processes. The code has been benchmarked against other hydrocodes 
(3§). A half-space target mesh was divided into two horizontal layers, the crust and the mantle, 
and cylindrical symmetry is assumed. The horizontal and vertical cell size was varied between 
0.5 x 0.5 and 1.5 x 1.5 km, depending on the impactor size. The choice of the cell size did 
not affect the results, as (a factor of 3) lower-resolution simulations (used for testing and not 
included here) also gave similar crust-mantle profiles and, most importantly, the same crustal 
thinning diameters. 

The pre-impact thickness of crustal materials was fixed at either 30 km (corresponding to 
the nearside) or 60 km (corresponding to the maximum thickness of the farside highlands). The 
impact speed was taken to be 10 or 17 km s -1 to accommodate a wide distribution of expected 
impact speeds during the epoch of basin formation (40), and the impactors (assumed to be made 
of dunite) ranged from 15 to 90 km in diameter in order to generate basins comparable in size 
to the observed lunar basins. These impact speeds could also represent moderately oblique inci- 
dence angles, because the vertical component of velocity vector mainly controls the final basin 
diameter (41). A lunar surface gravitational acceleration of 1.62 m s~ 2 was used. To allow 
faster computation, low-density material (<50-300 kg m 3 ) that largely represents vaporized 
material was removed from simulations, as it does not contribute to the basin formation pro- 
cess (7). The simulations were stopped 2 hours after impact. Temperatures beneath the basin 
were still elevated at that time, but the depth to the crust-mantle interface and the total crustal 
thickness are not expected to change markedly during further cooling (42). 

The material models for the crust, mantle, and impactor use the ANEOS-derived equation 
of state (EOS) tables, as well as strength (22), failure, and thermal softening models for basalt 
and dunite that were used in previous studies (24, 25) for modelling basin formation on the 
Moon and Mars (Table S2). These material parameters are similar to the ones used in other 
previous impact modelling studies that involved lunar mantle (42, 43). It was found in this 
study (by replacing basalt with granite) and other studies (42) that the choice of EOS for the 
crustal material does not make a substantial difference in the outcome of a basin-forming event. 

In this study, the role of temperature is restricted to its effect on the shear strength. The 
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material shear strength model in iSALE includes a description of thermal softening (22, 44): 

Y = Y c tanh [e(T m /T - 1)], (1) 


where Y c is is the cold shear strength, e is a material constant, T is temperature, and T m is the 
melting temperature. In this model, the strength decreases with temperature and is zero at, and 
above, the solidus. Because the solidus is not defined by the equation of state tables, we use the 
Simon approximation to fit the pressure dependence for the melting temperature of anhydrus 
basalt and KTB peridotite {45, 46) for the crust and mantle, respectively: 

T m = T mfi (p/ a + l) 1/c , (2) 

where T m 0 is the melt temperature at zero pressure, p is pressure, and a and c are material 
constants. 

Even though ANEOS is the best available representation of a multi-phase geological ma- 
terial, one limitation of the current ANEOS code for both materials is that it does not include 
the latent heat of melting at the transition between solid and liquid. For this reason, the temper- 
atures in the table above the solidus are overestimated when the rock is molten. However, as 
accurate melt volume calculations are not the focus of this study, this limitation is not important. 

Simulations employ the block-oscillation model of acoustic fluidization (46-48) to facil- 
itate crater collapse, which is important for cooler targets. A range of acoustic fluidization 
parameters were tested and varied until a basin forming in a cooler target collapsed into a final 
basin morphology with an acceptable basin depth (5-10 km) after the simulation ended. These 
parameters are also similar to the ACFL parameters employed in the simulations of Chicxu- 
lub crater collapse (49). Subsequent long-term cooling and relaxation of a basin over millions 
of years could cause the uplift of the complete basin for another few kilometers, essentially 
compensating for this depth (42). 

Although GRAIL has revealed that the lunar crust has an average porosity of 7-12%, for 
simplicity we did not employ iSALE’s porosity compaction model (23) in our simulations be- 
cause we do not expect porosity compaction to play an important role in large basin formation 
(as this is driven by uplift of deep mantle rocks with little porosity) and, as far as the data show, 
no major nearside-farside difference in porosity is observed. Including crustal pore-space com- 
paction would enhance shock heating of the crustal rocks in the basin center and dissipate the 
shock wave in the crust, possibly resulting in a slightly smaller crater. 

The input parameters to our simulations are provided in Table S2. The material model uses 
customized vertical thermal profiles of crust and upper mantle, as described in the following 
section. 
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Table S2. Target material parameters for iSALE-2D numerical models of impact basin formation. 


Description 

Crust 

Mantle and Impactor 

Equation of state 

Basalt ANEOS 

Dunite ANEOS 

Melt temperature at zero pressure (K), 7 m ( | 

1360 

1436 

Constant in thermal softening law, e 

0.7 

2.0 

Constant in Simon approximation (GPa), a 

4.5 

1.4 

Exponent in Simon approximation, c 

3.0 

5.0 

Poisson’s ratio 

0.25 

0.25 

Cohesion (damaged) (MPa) 

0.01 

0.01 

Coefficient of internal friction for material (damaged) 

0.6 

0.6 

Limiting strength at high pressure (damaged) (GPa) 

2.5 

3.5 

Cohesion (intact) (MPa) 

20 

50 

Coefficient of internal friction for material (intact) 

1.4 

1.5 

Limiting strength at high pressure (intact) (GPa) 

2.5 

3.5 

Minimum failure strain for low pressure states 

10" 4 

10” 4 

Increase in failure strain with pressure 

10- 11 

10" 11 

Pressure above which failure is always compressional (MPa) 

0.3 

0.3 

7 ^ constant 

0.004 

0.004 

7^3 constant 

230 

230 

Peak vibrational velocity as a fraction of the peak particle velocity 

0.1 

0.1 


3 Lunar thermal evolution modelling 

The three-dimentional (3D) spherical, finite-volume thermo-chemical convection code GAIA 
(2(5) was used to simulate the thermal evolution of the Moon. Our simulations, described by 
(20), follow previous work (15) by including an asymmetric distribution of heat sources as 
implied by orbital gamma-ray remote sensing measurements (50). To simulate the effects of 
the enhanced heat production from the Procellarum KREEP Terrane (13, 14) on the nearside 
hemisphere of the Moon, the equivalent of a 10-km-thick layer of KREEP basalt was placed 
in a spherical cap with an angular radius of 40°. For one model, the 10-km layer of enhanced 
heat production was placed just beneath a 40-km-thick crust (model PKT1), whereas for the 
other model, the equivalent of 10 km of KREEP basalt was distributed over the entire 40 km 
of crust (model PKT2). Our model has a bulk silicate abundance of 25 ppb uranium, which is 
comparable to the bulk silicate value on Earth. 

GAIA solves the standard hydrodynamics equations for an incompressible fluid, with an 
infinite Prandtl number under the Boussinesq approximation. Both core cooling and radioactive 
decay are included. The viscosity rj is taken to be Newtonian. Mantle melting is considered 
through the consumption of latent heat in the energy equation, and melt is assumed to leave the 
system instantaneously so that no latent heat of crystallization is released later. The solidus, T sol , 
and liquidus, Ti iq , are taken to be those of KLB-1 peridotite (51), and the density and depletion 
of the mantle (cumulative melt fraction) is tracked using a composition field. 

Two initial temperature profiles were considered for the lunar thermo-chemical evolution. 
For the first model, Ml, the initial temperature follows the solidus in the uppermost 350 km 
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of the mantle and an adiabatic gradient below (52). Conceptually, the region at the solidus 
represents the upper portion of the solidified magma ocean that did not convectively readjust, 
whereas the lower portion of the mantle underwent rapid convection. The second model, M2, 
used an adiabatic profile for the entire mantle (25), which corresponds to a well-mixed, initially 
convecting interior. In contrast to spherically symmetric thermal evolution models, for which 
the global cooling rate dictates the mantle behavior, the high heat production associated with the 
PKT dominates the thermal models here. The PKT heats the underlying mantle and gives rise to 
a shallow temperature peak for the nearside temperature profiles (Fig. S2). The farside is nearly 
unaffected by the PKT region, and its thermal evolution follows approximately a conductive 
state. On the nearside, regardless of the initial mantle temperature profile employed, models 
with the heat production enhancement below the crust (PKT1) give rise to a strong peak in 
temperature in the upper mantle. The model with the enhanced heat production uniformly 
distributed within the crust (PKT2) does not have such a peak as a result of the higher rate of 
heat loss to space, but it is still considerably hotter than the farside hemisphere. 
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Fig. S2. Vertical temperature profiles used in the impact simulations. Nearside and farside temperature 
profiles are shown with solid and dotted curves, respectively. Ml and M2 correspond to two different 
initial mantle temperature profiles; PKT1 and PKT2 include the enhanced heat production at the base 
of the nearside crust or distributed uniformly within the nearside crust, respectively. The temperature 
profiles for the upper mantle and crust on the farside are not affected by the nearside Procellarum KREEP 
Terrane. 
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The different temperature profiles used in our impact simulations correspond to the modeled 
thermal state after 0.5 billion years of evolution (i.e., 4 Ga). We also investigated an additional 
temperature profile for M2/PKT1 corresponding to 3.5 Ga, which is about 200 million years 
younger than Orientale. This additional model was included to show that the temperature effects 
persisted over a long time period (comparable to the duration of the basin-formation epoch) and 
how small changes in temperature profiles over this time would affect basin formation. 


4 Impact modelling results 

The parameters used in the iSALE impact simulations presented in Fig. 4 are shown in Tables S3 
and S4. The relation in Fig. 4 (reprinted with additional annotations and data in Fig. S3) between 
thin and thick crust for the same temperature profile is 

Ahin = 1.317 Awck, (3) 

where D is in km. This relation demonstrates that crustal thickness has little effect on the 
size of the final region of crustal thinning, as this relationship is very close to D thin / D thick ~1. 
Nevertheless, the absolute thickness of crustal material in the basin center, as well as the amount 
of excavated mantle material, depends on the pre-impact crustal thickness. 

The relations between crustal thinning on a hot nearside and that on the cool farside for 
different temperature profiles are: 


-Dfarside 

— 4 797 n 0 - 836 

^nearside) 

M1/PKT1 (4.0 Ga) 

(4) 

-^farside 

— 1 247 n 1 - 015 , 

^nearside) 

M2/PKT1 (4.0 Ga) 

(5) 

-^farside 

= 4.630 

M2/PKT1 (3.5 Ga) 

(6) 

-^farside 

= 0.796 

M1/PKT2 (4.0 Ga). 

(7) 


The differences in crustal thinning between the nearside and farside are greatest for the 
hottest nearside temperature profile at 4 Ga (M1/PKT1, shown in orange in Fig. S3). This rela- 
tion is very similar to that for the temperature profiles M2/PKT1 at 3.5 Ga (yellow in Fig. S3). 
This similarity demonstrates that our results are not critically dependent on the assumed time of 
basin formation, because even for a step as big as 0.5 Gyr the temperature effects persisted. In 
addition, Table S4 and Fig. S3 include impact modelling results for two different impact speeds, 
10 and 17 km s -1 . Under such a change in impact speeds, the relation between nearside and 
farside thinning diameters does not change, indicating that Eqs. (4)-(7) are largely independent 
of impact conditions. 
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Table S3. Impact simulation parameters used in Figs. 4 and S3 for the 
dependence on crustal thickness. L is the projectile size, v is the impact 
velocity, and D is the diameter of crustal thinning for a thick (60 km) or 
thin (30 km) crust. All thermal profiles correspond to the farside at 4 Ga. 
On the farside, M1/PKT1 is equal to M1/PKT2 and M2/PKT1 is equal 
to M2/PKT2, which means that the temperature increase characteristic 
for the nearside has no influence on the farside. 


L (km) 

v (km s 1 ) 

Temperature profile 

Ahick (km) 

Ahin (km) 

30 

17 

M1/PKT1 

155 

166 

30 

17 

M2/PKT1 

159 

167 

45 

17 

M1/PKT1 

244 

256 

45 

17 

M2/PKT1 

262 

265 

60 

17 

M1/PKT1 

334 

346 


Table S4. Impact simulation parameters used in Figs. 4 and S3 for the 


dependence on temperature. L is the projectile size, v is the impact 
velocity, and D is the diameter of crustal thinning on a thin nearside 
and thick farside crust. Temperature profiles are noted in time before 
present. 

L, km 

v, km s 1 

Temperature profile Gfarside, km 

-^nearside? km 

30 

17 

M1/PKT1 (4.0 Ga) 155 

315 

45 

17 

244 

433 

60 

17 

334 

572 

90 

17 

470 

797 

45 

17 

M2/PKT1 (4.0 Ga) 262 

359 

60 

17 

334 

450 

90 

17 

442 

610 

15 

17 

M2/PKT1 (3.5 Ga) 77 

154 

60 

17 

307 

538 

90 

17 

421 

684 

15 

10 

43 

110 

45 

10 

160 

328 

60 

10 

226 

420 

30 

17 

M1/PKT2, (4.0 Ga) 155 

184 

45 

17 

244 

312 

60 

17 

334 

420 
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Fig. S3. Dependence of impact basin size on target properties (from Fig. 4 with additional annotations 
and data). The diameter of crustal thinning for impact on a thin and hot nearside crust is plotted as a 
function of the corresponding diameter on a thick and cool far side crust. Simulations using the same 
temperature profile for the two hemispheres are shown in grey, whereas simulations with a hot nearside 
temperature profile resulting from the Procellarum KREEP Terrane are shown in other colors. Points of 
the same color correspond to simulations with different projectile diameters. The employed temperature 
profiles are given in Fig. S2. The dashed gray line marks a ratio of 1:1. 

The three vertical columns A, B, and C in Fig. S4 show series of snapshots from iSALE-2D 
impact simulations for the first 2 h after impact. All impacts were modelled with a 45-km- 
impactor at 17 km s -1 vertical impact velocity into: (A) 60-km-thick and cold farside crust, 
(B) 30-km-thick and cold farside, and (C) 30-km-thick and hot nearside crust. The temperature 
profile for A and B was M1/PKT2 (which is the same as M1/PKT1, as the PKT anomaly does 
not affect the thermal state of the farside) and for C it was M1/PKT2 (Fig. S2). Crust and 
mantle are represented by material concentration within the numerical mesh, which remains the 
same after the basin has cooled. The abscissa shows basin radius, in km, and the ordinate shows 
basin depth, in km, starting at the pre-impact surface taken as the zero-level. The values for the 
diameters of crustal thinning, D, for cases A, B, and C are in Tables S3 and S4. The first row 
of snapshots shows cratering initiation, the second row shows opening of the transient cavity 
and the latter is followed by a rapid crater floor rebound (third row). The fourth row shows the 
approximate maximum crater floor rebound, which is followed by collapse into a basin inner 
ring, as shown in the subsequent rows. The diameter of crustal thinning is marginally affected 
by different crustal thicknesses, considering that the diameter is measured as the radial distance 
from basin center at which the crustal thickness is first equal to the pre-impact crustal thickness. 
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Fig. S4. Snapshots from three iSALE-2D impact simulations, for the same impact conditions and differ- 
ent target properties. (A) Impact into thick crust corresponding to the cold farside highlands; (B) Impact 
into thin crust corresponding to the cold farside (e.g. South Pole-Aitken basin floor); (C) Impact into thin 
crust corresponding to the hot nearside hemisphere (PKT and mare-affected regions). 









This study was conducted only for a vertical impact scenario. Statistically more probable 
oblique impact scenarios would require 3D simulations that are computationally challenging. 
Because this work is a comparative study of basin formation on the nearside and farside hemi- 
spheres, and no variation of impact angle is expected between the hemispheres, we did not 
consider impact angle as an important parameter. However, we have performed a few oblique 
impact simulations in iSALE-3D (53-55) for both nearside and farside targets at 30° and 45° 
impact angles (measured from the horizontal). Although lowering the impact angle does sys- 
tematically reduce the size of the transient crater (as expected from crater scaling presented 
in (54)), and hence the amount of uplift and the diameter of the region of crustal thinning, this 
trend was observed for basins on both hemispheres and so does not affect our conclusions. 


5 Rescaling of the largest lunar impact basins 

With the relations in Eqs. (4)-(7), we determined the size that the nearside basins would have if 
they had formed on a cold farside crust (Table S5). This procedure is necessarily approximate as 
the subsurface temperature distribution is uncertain, because it depends on time, distance from 
the Procellarum KREEP Terrane, and the conditions that followed magma ocean crystallization. 
Regardless, we expect that our temperature profiles should compensate for the first-order differ- 
ences between the nearside and farside hemispheres. Considering that the Imbrium (341° E, 37° 
N), Serenitatis (19° E, 26° N), and C25 (351° E, 1 1° N) basins lie largely within the PKT, those 
basins were corrected using the hottest considered nearside temperature profile (M1/PKT1). 
The remainder of the nearside basins, which are all adjacent to mare deposits, were corrected 
using the coolest of the three considered nearside temperature profiles (M1/PKT2). 


Table S5. Diameter of crustal thinning D a hs. f° r the nearside and farside basins. Also listed for the near- 
side basins is the corresponding diameter D* if the basin had formed in crust with the same temperature 
profile as the farside. Entries in bold correspond to those shown in Fig. 2. 


Nearside basin 

^obs 


D * 


Farside basin 

-^obs 

M1/PKT1 

M2/PKT1 

M1/PKT2 

Imbrium 

680 

402 

497 

514 

Orientale 

418 

Serenitiatis 

607 

350 

444 

463 

Moscoviense 

319 

Crisium 

476 

262 

350 

370 

Freundlich-Sharanov 

308 

Smythii 

465 

255 

342 

362 

Coulomb-Sarton 

305 

Nectaris 

434 

235 

319 

339 

Fitzgerald-Jackson 

303 

Humorum 

380 

200 

280 

300 

Mendel-Rydberg 

286 

Asperitatis 

370 

194 

273 

293 

T22 

281 

Humboldtianum 

325 

166 

240 

260 

Apollo 

276 

Cruger-Sirsalis 

270 

133 

200 

219 

Hertzspmng 

252 

Humboldt 

235 

113 

174 

192 

Compton-Belkovitch 

230 

S chiller-Zucchius 

209 

98 

155 

173 

Korolev 

223 

C25 

235 

113 

174 

192 

Dirichlet-Jackson 

214 
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